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Abstract 

We present preliminary results from a search for the rare, flavor-changing neutral current decays 
B — >• Ki^i~ and B — > K*£~^£~, where i~^£~ is either an e~^e~ or fi^fj.~ pair. The data sample 
comprises (84.4 ± 0.9) x 10^ T{AS) BB decays (77.8 fb"^) collected with the BABAR detector 
at the PEP-II B Factory. For B — > K£^£~, we observe a signal with estimated significance of 
4.4cj and obtain B{B K£+£-) = (0.78+[5:^^;[5:J^) x 10^^ (averaged over £ = e and /i). For 
B K*£~^£~, we observe an excess of events over background with estimated significance of 



2.80-. We obtain B{B K*£+£-) = (LGStoH ± 0.28) x IQ-^ and the 90% C.L. upper limit 
B{B K*£+£~) < 3.0 X 10"^. 
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1 Introduction 



The flavor-changing neutral current decays B K£~^£~ and B K* {892)i~^ £~ , where is a 
charged lepton, are highly suppressed in the Standard Model, with branching fractions predicted 
to be of order 10~^ — 10~® |2|, The dominant contributions arise at one-loop level and include 
both electroweak penguin and box diagrams. Besides probing Standard Model loop effects, these 
rare decays are important because their rates and kinematic distributions are sensitive to new, 
heavy particles — such as those predicted by supersymmetric models — that can appear virtually in 
the loops i- 

Searches for B A'(*)£+£- decays have been performed by BABAR |], Belle §, CDF g, 
and CLEO 0]. Using a sample of 29 fb~^, the Belle Collaboration has observed a signal for 
the decay B K£+£- with branching fraction B{B Ke+i') = (0.75l|]|f ± 0.09) x lO'^^, 
averaged over electron and muon modes. The recently published BABAR result was based on 
a 20.7 fb-i sample and yielded the 90% C.L. upper limits B{B Kl+l') < 0.51 x 10"^ and 
B{B K*i~^£~) < 3.1 X 10~^. In the present study, we update our result with a total sample that 
is nearly four times larger than that used in the original BABAR analysis. 

We search for the following decays: B+ ^ K+£+i-, B° K'^l+i-, B+ K*+£+i-, and 
B° ^ K*°e+£-, where K*'^ K+tt', K*+ K^t^^ , K'^ tt+tt', and £ is either an e or /i. 
Throughout this paper, charge-conjugate modes are implied. 

2 Theoretical predictions 

The Standard Model processes for B K^*^i~^i~ include three contributions: the electromagnetic 
(EM) and Z penguin diagrams, and the W~^W~ box diagrams (Fig. [l|). Evidence for the EM 
penguin amplitude was first obtained by the CLEO experiment from the observation of the exclusive 
decay B K*j and later from the inclusive process B — > Xg'j, where Xg is any hadronic system 
with strangeness Since then, several experiments have observed these decays, and the Particle 

Data Group calculates the world averages B{B° K*^-i) = (4.3 ± 0.4) x 10"^ and B{B 
Xs-i) = (3.3 ±0.4) X 10"! 

The Standard Model rates for B K^i+l- and B Xs£^i~ are expected to be 25-80 
times smaller than those for B K*^ and B —>■ Xsj. For example, Ali et al. predict Q B{B — > 




w 

q q q q 

Figure 1: Standard Model diagrams for the decays B K^*k+i-. 
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XsC+e-) = (6.89±1.01)xl0-*^ and 13{B X.^u+zi") = (4.15±0.70) x lO'^ for the inclusive branch- 
ing fractions. Exclusive final states have larger theoretical uncertainties but are experimentally more 
accessible. AH et al. || predict B{B Ki+l') = (0.35 ± 0.12) x 10"^ for both e+e" and fi+fi' 
final states, B{B K*e+e~) = (1.58 ±0.49) x 10-^ and B{B K*fi+n~) = (1.19 ±0.39) x 10^^. 

Figure ^ shows the predicted rates as a function of q'^ = rn'j+^- for a number of theoretical 
models. The B — > K£~^£~ modes display a falloff in the rate from low values of q'^ (fast recoil 
kaon) to high values (slow recoil kaon). (Although B — s- Kj is forbidden by angular momentum 
conservation, B — > Ki~^i~ is allowed since q^ > 0.) The q^ distributions for B — > K*£^i~ behave 
quite differently, with large enhancements at low g^, but otherwise with the maximum rate at some 
intermediate value. The low q'^ enhancements are due to the EM penguin amplitude, which has 
a pole at q^ = 0, strongly enhancing B — > K*e~^e~ and, to a lesser extent, B — > K*fi~^n~. As a 
consequence, most Standard-Model-based predictions give a higher rate for B K*e~^e~ than for 
B — > K*fi~^n~. At values of q'^ above this pole region, the contributions from the Z penguin and 
the box diagram become very important. 

While the B K^*k+i- decays are rarer than the pure EM penguin processes, they have 
sensitivity to three Wilson coefficients C7, Cg, and Cio in the Operator Product Expansion Q, 
whereas the 6 — > 57 rate is sensitive mainly to the magnitude of a single Wilson coefficient, Cj. 
New physics could, for example, change the sign of C7 in S — > K*i~^i~; the resulting modification 
of the interference terms with the other amplitudes could enhance the rate by a factor of two over 
the Standard Model prediction. 

Measurements of the kinematic distributions in i? — > K*i~^i~ decays will eventually provide a 
probe of new physics with less theoretical uncertainty than the decay rates. Certain features of the 
lepton angular distribution as a function of can be predicted reliably and can be dramatically 
modified by new physics. In particular, the lepton forward-backward asymmetry in the dilepton 
rest frame is predicted to become zero at a particular value of q^ in the Standard Model, but the 
sign of the asymmetry and the location of the zero can be modified by new physics. The study of 
these angular distributions will require very large data samples and will be investigated by future 
experiments, including those at the LHC. 



3 The BaBAR detector and data sample 

The data used in the analysis were collected with the BABAR detector at the PEP-II asymmetric 
energy e^e~ storage ring at the Stanford Linear Accelerator Center. We analyze the data taken 
in the 1999-2002 runs, consisting of a 77.8 fb~^ sample accumulated on the T{AS) resonance, as 
well as 9.6 fb~^ taken at a center-of-mass energy 40 MeV below the T{4S) resonance peak to 
obtain a pure continuum sample. Continuum events include non-resonant e~^e~ qq production, 
where q = u, d, s, or c, as well as two-photon and t~^t~ events. The on-resonance sample contains 
(84.4 ± 0.9) X lO*' T{4:S) BB events. 

The BABAR detector is described in detail elsewhere [11 1. Of particular importance for this 



analysis are the charged-particle tracking system and the detectors used for particle identification. 
At radii between about 3 cm and 14 cm, charged tracks are measured in a five-layer silicon vertex 
tracker (SVT). Tracking beyond the SVT is provided by the drift chamber, which extends in radius 
from 23.6 to 80.9 cm and provides up to 40 track-measurement points. Just outside the drift 
chamber is the DIRC, a Cherenkov ring- imaging particle identification system. The DIRC provides 
charged particle velocity information through the measurement of the Cherenkov light cone opening 
angle. Cherenkov light is produced by charged particles as they pass through an array of 144 five- 
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Figure 2: Model predictions [|| for the decay rate as a function of cp' for (a) B — > Ke^e~ (b) 
B Kfi~^fj,~ (c) B — > K*e~^e~ and (d) — > K*fi~^fj,^. The solid histogram is based on the 
light cone QCD sum rules model of Ali et al., the dashed histogram on the QCD sum rules model 
of Colangelo et al., and the dot-dashed histogram on the quark model of Melikhov, Nikitin, and 
Simula. The distributions are normalized to the same area, so that any differences in the overall 
rates are not shown. 



meter- long fused silica quartz bars. The Cherenkov light is transmitted to the instrumented end of 
the bars by total internal reflection, preserving the information on the angle of the light emission 
with respect to the track direction. The DIRC is used for kaon identification in this analysis 
and is essential to our background rejection. Electrons are identified using an electromagnetic 
calorimeter comprising 6580 thallium-doped Csl crystals. These systems are mounted inside a 1.5 T 
solenoidal superconducting magnet to provide momentum measurement in the charged particle 
tracking systems. Muons are identified using the Instrumented Flux Return (IFR), in which resistive 
plate chambers (RPCs) are interleaved with the iron plates of the magnet flux return. 
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4 Analysis overview 



The most obvious analysis challenge is that in each of the eight signal channels there are back- 
ground processes with identical final-state particles that proceed through charmonium resonances. 
For example, the decay B K~^l~^£~ has backgrounds from B — > J/tp[-^ i~^i~)K~^ and B — >■ 
V'(25)(— > i~^i~)K'^ that must be carefully understood and vetoed. These backgrounds are par- 
ticularly difficult in the electron channels, where bremsstrahlung in detector material can easily 
cause the dilepton mass to fall well below that of the mass of the charmonium resonance. On the 
other hand, these backgrounds provide nearly ideal control samples for each channel, allowing us 
to make data vs. Monte Carlo comparisons and checks of selection criteria efficiencies in a way 
that carries over almost directly to the signal processes. Before searching for a signal in any newly 
acquired data sample, we validate our analysis by comprehensively studying the yields and shapes 
of distributions in the charmonium control samples. Although these processes have a narrow dis- 
tribution in q'^, they still allow us to confirm that the lepton identification efficiencies used in the 
Monte Carlo simulations are valid over a broad range of momenta. We have therefore studied the 
charmonium control samples in detail, not only to design our vetoes, but also to establish signal 
shape parameters and to verify Monte Carlo predictions of signal efficiencies. 

The most important non-charmonium background contributions are from BB backgrounds with 
either two real leptons or one real lepton and one hadron misidentified as a lepton (usually as a 
muon); continuum processes, especially cc events with a pair of D — > K^*^i~^i' decays or events with 
hadrons faking leptons; and small, but potentially very important contributions from a number of 
B decay modes with similar topology to the signal, such as B~^ — > D^vr"'' with K~tt~^, in 

which hadrons misidentified as leptons can create a false signal. Another example of a signal-like 
background is — > K*^ followed by conversion of the photon into e~^e~ in the detector material. 
This process is distinct from the decay B — > K*e~^e~ at low q^, where the photon undergoes internal 
conversion. While no single such background is large, it is important to establish that the sum of 
many processes with very small contributions does not produce a significant number of peaking 
background events. 

In the B Ki~^£~ modes, the dominant background is from continuum events. Much of this 
background can be suppressed using event shape variables that discriminate between the jet-like 
topology of continuum events and the much more spherical topology of signal (and other BB 
events). In i? ^ K*£'^l~ channels, the BB background is substantially larger than that from the 
continuum. Much of the BB background is due to pairs of primary semileptonic decays from the 
B and B mesons. Such events are characterized by large amounts of missing energy associated 
with the two neutrinos. We have constructed two main analysis variables, one for BB rejection 
{B likelihood) and one for continuum-background rejection (Fisher discriminant). Each of these 
variables combines information from several kinematic quantities. 

To prevent bias in the analysis, we optimize the event-selection criteria using (1) Monte Carlo 
samples both signal and background, (2) the off-resonance data, and (3) a "large sideband" region 
in the data that is outside of the fit region (defined below). Signal efficiencies were determined 
using the model of Ali et al. ^ , although a set of other models [||] was used to estimate a systematic 
uncertainty on the efficiencies due to model dependence. None of our models for the signal efficiency 
include small effects arising from interference between the electroweak penguin amplitudes and the 
amplitudes for the decays involving charmonium intermediate states. The interference is construc- 
tive on the lower side of the J and ip{2S) charmonium resonance masses but is destructive on 
the upper side. The effect of this interference on our yields is also suppressed by the charmonium 
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vetoes. 

After event selection, we perform a two-dimensional unbinned maximimum likelihood fit to the 
joint distribution in 



mES = jK'-iT.P*) 



i 

where is the beam energy in the e^e~ rest (CM) frame, p* is the CM momentum of daughter 
particle i of the B meson candidate, and mj is the mass hypothesis for particle i. For signal events, 
rriEs peaks at the B meson mass with a resolution of about 2.5 MeV/c^ and /S.E peaks near zero, 
indicating that the candidate system of particles has total energy consistent with the beam energy 
in the CM frame. 

Figure ^ shows three regions that we have defined in the /S.E vs. ttt-es plane. The events are 



taken from a GEANT4 [12| simulation of the process B^ K^e^e~ . The small box where most 
of the events lie is the nominal signal region. This region is used to indicate where most of the 
signal is expected, but signal events are extracted over the full fit region, which is much larger. The 
fit region is concealed during the process of determining the event selection criteria. The region 
outside the fit region is called the large sideband. It is not concealed but is instead used to perform 
data vs. Monte Carlo comparisons in a region that is dominated by combinatorial background. This 
region, together with samples of Ke^ and K*eiJL events, provides a way for us to study combinatorial 
backgrounds without the possibility of biasing the sample of events in the fit region. In the fits, the 
background normalization is allowed to float, as are other parameters describing the background 
shape. Thus, while we use Monte Carlo samples to motivate our background parametrizations, the 
parameter values themselves are not fixed from the Monte Carlo simulation. 



5 Event selection 

We select events that have at least four charged tracks, the ratio R2 of the second and zeroth Fox- 



Wolfram moments [13| (calculated using charged tracks and unmatched calorimeter clusters) less 
than 0.5, and two oppositely charged leptons with labratory-frame momenta p > 0.5 (1.0) GeV/c 
for e (^) candidates. The leptons are identified with strict criteria to minimize the misidentification 
rates for hadrons. The typical probability for a pion to be misidentified as an electron is 0.2% to 
0.3%, depending on momentum; the probability for a pion to be misidentified as a muon is 1.3% 
to 2.7%. Electron-positron pairs consistent with photon conversions in the detector material are 
vetoed. In the B Ke^e~ modes, e^e~ pairs at low are vetoed independent of the distance of 
the vertex from the beam axis, but ui B ^ K*e'^e~, where we expect substantial rate at low q^, 
we veto only conversion candidates that occur at a radius of at least 2 cm, where detector materials 
are present. We require charged kaon candidates to be identified as kaons with strict criteria and 
the charged pion in K* — > Kir not to be identified as a kaon. The use of strict criteria for kaon and 
lepton identification is necessary to suppress peaking backgrounds. For B K*i^i~ , we require 
the mass of the K* candidate to be within 75 MeV/c^ of the mean K*(892) mass. Kg candidates 
are reconstructed from two oppositely charged tracks that form a good vertex displaced from the 
primary vertex by at least 1 mm. The mass of the 'k~^tt~ system must be within ±9.3 MeV/c^ of 
the nominal Kg mass. 
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Figure 3: The distribution of Monte-Carlo-simulated B'^ K'^e^e~ events in the rriEs vs. AE 
plane. The nominal signal region, fit region, and large-sideband region are shown. The full fit 
region in the data is concealed ("blinded") until the event selection procedure is defined and the 
data validation steps are performed. The nominal signal region is used only for the optimization 
of the event selection criteria using Monte Carlo samples; the extraction of the signal yields from 
the data is based on an analysis of the entire fit region. 



The decays B ^ J/ipi-^ i+i-)K'-*^ and B ip{2S){-^ i+i-)K'^*^ have identical topologies to 
signal events. These backgrounds are suppressed by applying a veto in the AE vs. m£+£- plane 
(Fig. 1^. This veto removes charmonium events not only with reconstructed m^+g- values near the 
nominal charmonium masses, but also events that lie further away in m£+£- due to photon radiation 
(more pronounced in electron channels) or track mismeasurement. Removing these events — not only 
from the signal region, but also from the full fit region — simplifies the description of the background 
shape. Charmonium events can, however, escape this veto if one of the leptons (typically a muon) 
and the kaon are misidentified as each other. If reassignment of particle types results in a dilepton 
mass consistent with the J/ip or ip{2S) mass, the candidate is vetoed. There is also significant feed- 
up from B J/ijjK and B — > ij){2S)K into B — > K*i~^i~ , since energy lost due to bremsstrahlung 
in i? — > J/ipK can be compensated for by including a random pion. If the Ki'^i~ system in a 
B K*i~^i~ candidate is kinematically consistent with B J £^i~^)K, assuming that 
the photon (which is not directly observed) was radiated along the direction of either lepton, then 
the candidate is vetoed. In modes with muons we also apply a veto against backgrounds from 
B Dtt, with D K'-*\. We reassign the particle hypothesis for the muon candidate with 
charge consistent to form a D when combined with the K^*'^ to a pion and veto the candidate if 
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Figure 4: Charmonium veto in the AE" vs. m£+£- plane for (a) B K^*^e~^e and (b) B — > 
. Hatched regions are vetoed. The dots correspond to a Monte Carlo simulation of 
B J/ipi-^ i+i-)K and B V(25')(^ i+l-)K. Most signal events lie in the AE region 
between the horizontal lines. 

the mass of the muon-X^*) system is consistent with the D mass. 

Apart from the charmonium vetoes, we analyze the full range of g^. In the B — > Ke~^e~, 
B — > Kii'^fi~, and B K*e~^e~ modes, we have good efficiency over most of the range, while in 
the B — > K* ^'^ n~ mode the efficiency is highest at intermediate and high values. This efficiency 
variation results from the interplay between the lepton identification efficiency and the reconstruc- 
tion efficiency of the recoil hadron as a function of momentum. The average lepton momentum 
increases as increases; since our minimum momentum requirement is lower for electrons than 
muons, the region of good lepton efficiency generally extends to lower values of q^ for electrons. 
However, the momentum of the recoil hadron decreases as q^ increases, which produces the opposite 
effect in the K^*"^ efficiency. 

Continuum background from non-resonant e^e~ — qq production is suppressed using a Fisher 
discriminant which is a linear combination of the input variables with optimized coefficients. 
The variables included are R2'-, cos 6b, the cosine of the angle between the B candidate and the 
beam axis in the CM frame; cos 6t, the cosine of the angle between the thrust axis of the candidate 
B meson daughter particles and that of the remaining particles in the CM frame; and mxi, the 
invariant mass of the X-lepton system, where the lepton is selected according to its charge relative 
to the strangeness of the K^*\ The variable rriKe helps discriminate against background from 
semileptonic D decays, for which mj^{,)^ < mo. The Fisher discriminant is optimized separately 
for each decay channel. 

Combinatorial background from BB events is suppressed using a signal-to- -B-B likelihood ratio 
that combines candidate B and dilepton vertex probabilities; the significance of the dilepton separa- 
tion along the beam direction; cos 9b', and the missing energy, -Bmissj of the event in the CM frame. 
The variable -Bmiss provides the strongest discrimination against BB background, since events with 
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semileptonic decays usually have significant unobserved energy due to neutrinos. We construct a 
likelihood variable because the shape of the B vertex probability distributions are not well suited 
to a Fisher discriminant. A separate B likelihood variable is constructed for each channel. 

For each final state, we select at most one combination of particles per event as a -B signal 
candidate. If multiple candidates occur in the fit region, we select the candidate with the greatest 
number of drift chamber and SVT hits on the charged tracks. In the B — > Ki~^l~ modes less than 
1% of the events have multiple candidates and in the B — > K*£'^i~ modes about 10% of the events 
have more than one candidate. This criterion is designed to be negligibly biased with respect to 
the kinematic quantities used in our fit. 



6 Control sample studies 

To check the ability of the Monte Carlo to correctly simulate the detector response, we study the 
charmonium decays B — > J/il)K^*^ and B — > iIj{2S)K^*\ whose decay distributions and branching 
fractions are reasonably well known. Figures ^, ^, and ^ compare the distributions of mi+i- , 
lepton energy in the T(45') CM frame, B likelihood variable, and Fisher continuum suppression 
variable for these signal- like control samples. In each case the distributions are absolutely normal- 
ized, so that we are comparing not only the agreement between shapes, but also the agreement 
between event yields. Table || summarizes the event yields in the nominal signal boxes for data 
and the Monte Carlo prediction in the control samples. We observe good agreement in the yield 
between data and Monte Carlo simulation in all channels. The agreement in the tails of the dilepton 
mass distributions indicates that the placement of our charmonium vetoes is based on a reliable 
simulation. The lepton-energy distributions are well simulated, indicating that the Monte Carlo 
efficiencies should be well simulated at other values of (f'. The selections on the B likelihood and 
Fisher continuum-suppression variables are rather loose, and the agreement between Monte Carlo 
simulation and data in the control samples indicates that the Monte Carlo efficiencies should be 
reliable. We discuss these issues quantitatively in the discussion of systematic errors (Sec. ^). 

The large sidebands provide a useful tool for monitoring our simulation of the combinatorial 
backgrounds. In contrast to the charmonium control samples, the events here are produced by 
a wide variety of processes, including both BB and continuum. We compare the distributions 
observed in data with the sum of BB Monte Carlo simulations and continuum samples obtained 
from off-resonance data. Although the off-resonance sample is statistically limited, it gives a correct 
picture of the complicated mix of qq^ t~^t~ , and two-photon events. Figure ^ shows the distribution 
of the B likelihood variable in the large sideband. Both the shapes and the normalizations are 
reasonably well predicted. There are large fluctuations in the off-resonance data sample, which is 
scaled by a factor of about 8 to correspond to the integrated luminosity of the on-resonance data 
sample. 

For each data taking period, we check the control samples in the manner described above. We 
then "unblind" the data in the fit region. The distributions of the data for each mode are shown in 



Fig. IC . In several of the modes, a small number of events is present in the nominal signal region. 



For the events in the nominal signal region we show the dilepton mass distribution in Fig. 11. The 
prediction based on our simulation of the signal is also shown. 
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Figure 5: Comparison of the m^+^- distributions between data and Monte Carlo simulation for the 
charmonium control samples. The points with error bars show the data, and the solid histograms 
show the prediction of the charmonium Monte Carlo simulation. All of the analysis selection criteria 
have been applied except for the charmonium veto, which is reversed. The events are also required 
to be in the nominal signal region in the toes vs. IS.E plane. The large tails in the e"'"e~ modes are 
due to photon radiation. 
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Figure 6: Comparison of the lepton-energy distributions between data and Monte Carlo simulation 
for the charmonium control samples. The points with error bars show the data, and the solid 
histograms show the prediction of the charmonium Monte Carlo simulation. All of the analysis 
selection criteria have been applied except for the charmonium veto, which is reversed. The events 
are also required to be in the nominal signal region in the mss vs. /S.E plane. 
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Figure 7: Comparison of the B likelihood distributions between data and Monte Carlo simulation 
for the charmonium control samples. The points with error bars show the data, and the solid 
histograms show the prediction of the charmonium Monte Carlo simulation. All of the analysis 
selection criteria have been applied except for that on the B likelihood and the charmonium veto, 
which is reversed. The events are also required to be in the nominal signal region in the rriES vs. AE 
plane. The vertical lines show the minimum B likelihood values allowed for the corresponding signal 
modes. 



18 




in 
o 



0) 

> 

LU 



-5 5 

Continuum Fisher 



10 




-5 5 

Continuum Fisher 



10 



50 
40 
30 
20 
10 




1 

- A 


K^e^e- _ 


7\ 





o 



0) 

> 

LU 



-10 



5 
Continuum Fisher 



10 



25 
20 
15 
10 
5 






1 




1 


1 


- / 


1 







-10 -5 5 

Continuum Fisher 



10 




o 



0) 

> 

LU 



-5 5 

Continuum Fisher 



10 




-5 5 

Continuum Fisher 



10 




-5 5 

Continuum Fisher 



5 5 

Continuum Fisher 



10 



Figure 8: Comparison of the Fisher continuum-suppression variable distributions between data 
and Monte Carlo simulation for the charmonium control samples. The points with error bars 
show the data, and the solid histograms show the prediction of the charmonium Monte Carlo 
simulation. All of the analysis selection criteria have been applied except for that on the Fisher 
and the charmonium veto, which is reversed. The events are also required to be in the nominal 
signal region in the mES vs. /S.E plane. The vertical lines show the maximum Fisher values allowed 
for the corresponding signal modes. 
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Figure 9: Large sideband comparison of the B likelihood variable distributions between data 
and BB Monte Carlo simulation plus scaled off-resonance data. The points with error bars show 
the on-resonance data, the shaded histograms show the scaled off-resonance data, and the open 
histograms show the sum of the scaled off-resonance data and the BB Monte Carlo sample. All of 
the analysis selection criteria have been applied except for that on the BB likelihood. The position 
of this cut is indicated by vertical lines; selected events are required to have B likelihood above 
this value. 
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Figure 10: Distribution of the data in the AE vs. 
rectangles show the nominal signal boxes, but the sij 



mEs plane (fit region) for each channel. The 
;nals are extracted from a two-dimensional fit. 
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Figure 11: The dilepton mass distribution for data (points with error bars) and signal Monte Carlo 
(histogram) for events in the nominal signal region. The data include some amount of background, 
so the distributions are not expected to agree perfectly. However, the data do not cluster around the 
charmonium veto regions; such clustering would be a possible indication of charmonium background 
events escaping the veto. 
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Table 1: Comparison of the J/-0 control sample event yields in data and Monte Carlo simulations 
for each of the modes. All selection criteria have been applied, except that the J/ijj veto is reversed, 
and the nominal signal region is selected in mES and /S.E. The Monte Carlo samples have been 
scaled to the equivalent number of T(4S') decays in the corresponding data sample. Errors on the 
ratios of data yields to yields from the Monte Carlo simulation are due to our sample statistics 
only, and do not incorporate branching fraction uncertainties. 



Mode (J/ip^l+r) 


MC Yield 


Data Yield 


Data/MC(%) 


K^e+e- 


911.8 


939 


103.0±3.8 


B^ K^fi+fi- 


669.1 


623 


93.1±4.1 


B^ K^e+e- 
B° K'^fi+fi- 


271.6 


265 


97.6±6.0 


195.0 


191 


98.0±7.1 




494.5 


523 


105.8±4.9 


BO K*^fi+i^- 


353.7 


412 


116.5±6.2 




158.6 


144 


90.8±7.7 


B± ^ K*^fi+fi- 


106.4 


111 


104.4±10.2 


Combined by hadron final state 






B^ K^i+e- 


1580.9 


1562 


98.8±2.8 


B^ K^£+£- 


466.5 


456 


97.7±4.6 


bo ^ K*°£+e- 


848.2 


935 


110.2±3.9 


B± ^ K*^£+i- 


264.9 


255 


96.3±6.2 


Combined by lepton 


final state 






e+e~ modes 


1836.4 


1871 


101.9±2.6 


H~^fi~ modes 


1324.1 


1337 


101.0±3.0 


all modes 


3160.5 


3208 


101.5±1.9 



7 Fits and yields 

We extract the signal and background yields in each channel using a two-dimensional extended 
unbinned maximum likelihood fit in the region 

defined by mES > 5.2 GeV/c^ and |A^j < 0.25 GeV 
("fit region"). The two-dimensional signal shapes, including the effects of radiation on the AE 
distribution and the correlation between tties and AE, are obtained by parameterizing a GEANT4 



Monte Carlo simulation of the signal using the product of two Crystal Ball shapes |15], one for 
rriES and one for AE. The Crystal Ball shape smoothly matches a Gaussian core with a tail that is 
used to describe radiative effects. For w-es we use a single Gaussian core, while in AE we describe 
the core using the sum of two Gaussians with separate mean and width. The parameters in the 
Crystal Ball shape for m-ES are allowed to have a quadratic dependence on AE, which adequately 
takes into account the correlation between mEs and AE. Small shifts in the mean values of rriES 
and AE are applied to the signal shapes based on studies of these quantities in the charmonium 
control samples. 

The combinatorial background is described by the function 



/(A£:,mEs) = Ne^^^m^sjl - '-^e\' / , (1) 
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where is a normalization factor, and s and are free parameters determined in the fit to the 
data. The part of this function that describes the mss distribution is a standard form known as 



the ARGUS function |16|. 



For systematic studies we also use a more general form for the background shape, where the 
ARGUS slope parameter, ^, is allowed to be a function of AE: 

C = aAE)=^o + ^iAE + ^U^Ef. (2) 

In each channel we also include a term to describe the estimated amount of background that 
peaks under the signal. Peaking backgrounds can arise from B K*^ followed by a photon 
conversion in the detector material or from decays such as B ^ Ktttt in which the pions are 
misidentified as leptons. Our estimates for these backgrounds are obtained from specific studies of 
individual processes and make use of hadron misidentification rates measured in data. The total 
estimated peaking background for each mode is listed in Table ^. 

Table 2: Summary of peaking backgrounds used in the fits. 



Mode 


















B^^ 

















Peaking background (events) 
~n} 



O.Olo.o 
0.7 ±0.7 

0.5 ±0.5 

1.2 ± 1.2 

1.0 ± 1.0 



As an example, we present the fits for the B — > K*^e^e^ mode. Figure ^ shows the ttt-es 
distributions in three slices of AE that span the fit region: AE < 0.11 GeV, —0.11 < AE < 0.05 
GeV, and AE > 0.05 GeV. We observe a peak in mss around the B mass for the range —0.11 < 



AE < 0.05 GeV, but not in the upper or lower AE sidebands. Figure |l^ shows the distributions 
of the data in AE for three slices of rriEs: "ies < 5.24 GeV/c^, 5.24 < tties < 5.273 GeV/c^, 
mES > 5.273 GeV/c^. 



Figure 14 shows the projections of the fits for each channel onto m^s for AE in the nominal 



signal regions. Figure 15 shows the projections of the fits for each channel onto AE for rriEs in 
the nominal signal regions. We observe small excesses of events over the backgrounds in several 
channels. Table |^ summarizes the results for each channel. 

We also perform combined fits in which the branching fractions in the four B — s- Ki^i~ channels 
are constrained to each other, as are the branching fractions for the four B — > K*i~^i~ channels 
(Fig. |l6|). Due to the larger pole at low values of mi+£- in B ^ K*e~^e~ than in B —> K* fi^ ji~ , 
the expected rates for these channels differ and are constrained to the ratio of branching fractions 
B{B K*e+e-)/B{B K* 11+ ^j.-) = 1.2 from the model of AH et al. 0. The extracted yield 
corresponds to the electron mode. The branching fractions from the combined fits, with statistical 
errors only, are 

B{B Kl+r) = (0.78l[5:i^(stat.)) x 10-^ 
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Table 3: Results from individual fits to B ^ K^*'^£~^i~ modes. The columns are, from left to 
right: (1) the fitted signal yield; (2) the contribution of the background to the error on the signal 
yield, expressed as an effective background yield; (3) the signal efficiency, e (not including the 
branching fractions for K*, K^, and K^, decays); (4) the multiplicative systematic error on the 
selection efficiency, {AB/B)e', (5) the additive systematic error from the fit, (A;B)fit; and (6) the 
branching fraction central value (B). 



Mode 



Signal 
Yield 



Effective e (AB/B), {AB){,t B 
Bkgnd. (%) (%) (/10-6) (/10-6) 



B+ 


K+e+e- 


14.4+^-^ 


2.2 


17.5 


±6.8 


+U.i4 
-0.21 


n 08+0.34+0.16 
^■^"-0.28-0.22 


B+ 






2.2 


9.2 


±6.6 


+0.09 
-0.08 


n nfi+0-30+0.09 

"-'•'-'"-0.17-0.08 


B° 


K°e+e- 


1 q+2.6 


1.5 


18.6 


±7.9 


±0.14 


n 9/1+0.49+0.14 
"J-^^-0.32-0.15 


B^ 


KV+M" 


3.6i:? 


0.9 


9.4 


±7.7 


±0.24 


1.33lo78 ±0-26 


B^ 




10.6+|-3 


5.5 


10.6 


±7.6 


+0.46 
-0.47 


1 7o+0:87+0.48 
'O-0.72-0.49 


B^ 




3.4tg 


4.6 


6.1 


±9.3 


±0.38 


0.99+J-^^ ± 0.39 


B+ 


K*+e+e- 


0.3+li 


6.3 


10.3 


±9.5 


+0.69 
-0.72 


n 1 t;+i:87+0.69 
"J- -'^^-1.16-0.72 


B+ 




3.6ta 


3.4 


5.2 


±11.1 


±1.80 


3Mtli\ ± 1.84 



B{B^K*i+r) = (1.68t^;^^(stat.)) X 10"^ 



8 Systematic uncertainties 

Systematic errors arise from uncertainties on the efficiency and on the number of BB pairs in the 
data sample (both of which are multiplicative errors) and on the yields extracted from the fits 
(additive systematic errors). 

Table ^ summarizes the multiplicative systematic uncertainties. The systematic uncertainty 
from the continuum Fisher and BB likelihood selection criteria are determined from the differences 
of efficiency for these criteria in data and Monte Carlo simulations for charmonium control samples 
(see Figures |^ and |8|). The uncertainties on the efficiencies due to model-dependence of form factors 
are taken to be the full range of variation obtained from different theoretical models [j^. The total 
multiplicative systematic uncertainty is the sum in quadrature of the individual sources, with the 
exception of the tracking efficiency uncertainties for leptons and hadrons, which are taken to be 
100% correlated. The total uncertainty ranges from 7-11%; the largest individual source is the 
theoretical model dependence of the signal efficiency, which ranges from 4-7%. The combined 
multiplicative systematic errors are also listed in Table ^. 

The additive systematic errors are uncertainties in the signal yield from the fit. These arise from 
three sources: (1) uncertainties in the signal shapes, (2) uncertainties in the background shapes, and 
(3) uncertainties in the amount of peaking backgrounds. Note that because the background shapes 
and yields float in the fit, much of the uncertainty associated with the background is automatically 
incorporated in the statistical uncertainty of the signal. A systematic uncertainty in the background 
shape is evaluated by (1) fixing the mss slopes to the values obtained from Monte Carlo simulations 
rather than allowing the slopes to float in the fit and (2) performing alternative fits in which the 
TTT-ES slope is allowed to have a quadratic dependence on AE, as described above. 
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m^^ (GeV/c ) 



Figure 12: Fit to the data for the mode K*^{^ K'^tt )e^e . The three histograms show 

the m-Es distributions in three shces of l^E: AE < -0.11 GeV, -0.11 < AE < 0.05 GeV, and 
> 0.05 GeV. 
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Figure 13: Fit to the data for the mode — > > K~^7r~)e~^e~ . The three histograms show 

the AE distributions in three slices of rriES- "lES < 5.24 GeV/c^, 5.24 < ttt-es < 5.273 GeV/c^, 
n^ES ^ 5.273 GeV/c^. The peak in AE near zero in the lower histogram is correlated with the 
peak in rn-ES around the B mass. 
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Figure 14: Projections of the individual fits onto the tties distributions for the shce —0.11 < Ai? < 
0.05 GeV (electron modes) or —0.07 < /S.E < 0.05 GeV (muon modes). 
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Figure 15: Projections of the individual fits onto the AE distributions for ttt-es in tlie nominal 
signal regions. 
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Figure 16: Summed distributions for the four B — > Ki~^i~ and the four B — > K*i~^i~ channels, 
showing the projections of the combined fit. In this combined fit, the yields in the B Ki~^£~ 
channels are constrained to each other, as are the yields in the B — > K*i~^£~ channels (see text). 
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Table 4: Multiplicative systematic uncertainties. 



Systematic 


A e e 




Kse e 




A e e 




i\ e e 




irk en. (e, fi) 


±1.6 


±1.6 


±1.6 


±1.6 


±1.6 


±1.6 


±1.6 


±1.6 


rjiectron lu 


-1-9 7 




-1-9 7 




-1-0 7 




-1-9 7 




Muon ID 




±2.0 




±2.0 




±2.0 




±2.0 


K,n ID 


±2.0 


±2.0 






±4.0 


±4.0 


±2.0 


±2.0 


Trkeff. (A.tt) 


±1.3 


±1.3 


±2.6 


±2.6 


±2.6 


±2.6 


±3.9 


±3.9 


K'^ eff. 






±3.2 


±3.2 






±3.2 


±3.2 


SB Counting 


±1.1 


±1.1 


±1.1 


±1.1 


±1.1 


±1.1 


±1.1 


±1.1 


Fisher 


±1.5 


±1.5 


±1.5 


±1.5 


±1.5 


±1.5 


±1.5 


±1.5 


BB likelihood 


±2.5 


±2.5 


±2.5 


±2.5 


±2.5 


±2.5 


±2.5 


±2.5 


Model dep. 


±4.0 


±4.0 


±4.0 


±4.0 


±4.0 


±7.0 


±4.0 


±7.0 


MC statistics 


±1.1 


±1.1 


±1.1 


±1.1 


±1.9 


±1.6 


±1.0 


±1.9 


Total 


±6.8 


±6.6 


±7.9 


±7.7 


±7.6 


±9.3 


±9.5 


±11.1 



The K~^e~^e~ channel has one additional systematic uncertainty included. In this channel 

we observe two events with the following property: if the electrons are each assigned pion masses, 
then the system is consistent with coming from the decay D^Tr~^; — > K~7r~^. The 

predicted number of such events, based on measured probabilities for pions to be misidentified 
as electrons, is only about 0.08 event. In spite of this small estimated probability we take into 
account the possibility that these events might in fact be background by including an asymmetric 
systematic error corresponding to the subtraction of two events. This issue does not arise in the 
muon channels, where we have an explicit veto on this background to handle the much higher 
probability for a pion to be misidentified as a muon. 

Table |5| summarizes the additive systematic errors, expressed directly in terms of errors on 
signal event yields. Table ^ gives the total additive systematic errors propagated to the branching 
fractions. For the combined fits for B — > K£^i~ and B K*£~^i~ the systematic errors are 
calculated by weighting the systematic errors from the individual modes based on the expected 
signal from the efficiency and branching fractions. 



Table 5: The sources of additive systematic uncertainty, expressed directly as uncertainties in the 
signal yield. 



Systematic 


K+e+e- 




Aiie+e- 












Signal shape 


±0.3 


±0.0 


±0.3 


±0.1 


±0.4 


±0.2 


±0.4 


±0.3 


Comb bkg. shape 


±1.8 


±0.1 


±0.7 


±0.3 


±2.6 


±0.3 


±1.3 


±1.4 


Peaking bkg. 


+0.0 
-2.2 


±0.7 


+0.0 
-0.1 


±0.5 


+0.2 

-0.5 


±1.2 


+0.1 
-0.4 


±1.0 



9 Discussion 

This analysis has yielded preliminary evidence for the decay B Ki'^£~. The signal events are 
primarily in the i?"^ — > K'^e'^e" channel, although there is also an excess in B^ Kgfi'^fi". 
The higher event yield in electron over muon channels is not unexpected, due to the substantially 
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higher detection efficiency for electrons than muons in the BABAR detector. The significance of this 
signal, based purely on the statistical errors, is 5. 60"; when systematic uncertainties are applied, 
the estimated significance becomes 4. 4a. We conclude that a signal for B Ki~^i~ has been 
observed. The measured branching fraction is consistent with most predictions, although the 
central value of B{B Ki+i') = (0.78+o;2otai8) x 10"^ is higher than the Ah et al. prediction 
of B{B K£~^£~) = (0.35 =b 0.12) x 10~^. The result is also consistent with Belle's measurement, 
B{B K£+£-) = (0.75l|5;^f ± 0.09) x 10"^. The result from this 77.8 fb^^ sample is also higher 
than the 90% C.L. upper limit obtained from our original analysis of 20.7 fb~^; the larger data set 
includes the smaller data set. We have made a conservative evaluation of the consistency in our 
results as follows: We evaluate the B Ki~^l~ branching fraction using only the data acquired 
after that used in the published result. Assuming a true value equal to the result obtained in this 
newer dataset, we find a 3% probability {1.9a) that we would have observed a yield smaller than 
our published result from the earlier dataset. 

We also observe an excess of events in — > K*i^i~ , with a significance including systematic 
uncertainty of 2.8a. We obtain a central value of B{B K*e+£-) = (1.68l[J|^ ± 0.28) x IQ-^ 
and a 90% C.L. upper limit B{B K*£~^i~) < 3.0 x 10"^. The central value is within the range 
expected in the Standard Model, although the uncertainties are still large. 

10 Conclusions 

We present preliminary results of our search for B K(*k+£- usmg a sample of 84.4 x 10^ BB 
pairs produced at the T{4S). We observe a signal for B K£~^£~ with 

B{B^K£^£-) = (0.78l0:2^;0:ii) X 10-6_ 
The significance is estimated to be A.4a. For B — > K*£'^£~ we obtain 

B{B K*£+£-) = (1.68+|]:| ± 0.28) x 10-^ 

using the constraint B{B K*e^e~)/B{B K*fi^fi~) = 1.2. For the combined B K*£^£~ 
result we quote the branching fraction corresponding to the electron channel. The significance of 
the B K* £'^ £~ signal is estimated to be 2.8cr and we quote a 90% C.L. limit for this mode of 

B{B K*£+£-) < 3.0 X 10"^ 
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